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Abstract 

A method of reducing the problem of the calculation of tree mul- 
tiparticle cross sections in </? 4 theory to the solution of a singular clas- 
sical Euclidean boundary value problem is introduced. The solutions 
are obtained numerically in terms of the decomposition in spherical 
harmonics, and the corresponding estimates of the tree cross sections 
at arbitrary energies are found. Numerical analysis agrees with ana- 
lytical results obtained earlier in the limiting cases of large and small 
energies. 



1 Introduction 

The main theoretical tool in quantum field theory is presently the pertur- 
bation theory. Perturbative calculations provided the majority of the exper- 
imentally checked results. Therefore, limits of applicability of perturbative 
calculations are of considerable interest. On the one hand, there exist pro- 
cesses related to complex vacuum structure in gauge theories and nontrivial 
classical solutions to field equations which cannot be described by perturba- 
tion theory. On the other hand, even in the topologically trivial sector at 
relatively low energies, processes which are described poorly by perturbation 
theory also exist. 
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Perturbative calculations provide reliable results only in weakly coupled 
models where the expansion parameter — dimensionless coupling constant — is 
much smaller than unity. But even in such theories situations are possible in 
which other competing small (or large) parameters exist. A typical example 
is a process with large number n of particles in the final state (n being of the 
order of the inverse coupling constant A -1 ). 

In conventional perturbation theory even above the topologically trivial 
vacuum, the naive estimate of the amplitude gives the factorial dependence 
n! on the multiplicity of the final state. This enhancement can in principle 
overcome the suppression due to the powers of the coupling constant. At 
the tree level, it is possible to find an exact expression for the amplitude 
of creation of n real particles by one virtual particle in the theory with 
lagrangian density 

c = \{d, v f - - \^ (i) 

(the mass is set equal to one) in special kinematics, namely, when all particles 
have zero spatial momenta 

n-l 

This result points towards complete breakdown of the usual perturbative 
calculations at n > A -1 because it contradicts unitarity of the theory. 

Thus, some non-perturbative method is required for the calculation of 
these cross sections. The limit we are interested in is 

A — ► , Xn = fixed , e = fixed , (3) 

where e = (E — n)/n is the average kinetic energy of the outgoing particles 
in the centre of mass frame. Existing perturbative calculations [|2], |3| strongly 
suggest that in this limit the total cross section has the exponential form, 

<Ti^ n ~ exp f ip(An,e)^ . (4) 

This form implies the semiclassical calculability of the cross sections. A 
method to obtain the exponent F(Xn,e) in all loops was formulated in ref. 
[[|], which required solution of a certain classical boundary problem in com- 
plex time. At small An one needs only to solve purely Euclidean equations 
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with special boundary conditions. In conventional perturbative approach this 
limit corresponds to the contribution of tree graphs, that gives the following 
dependence on A: 



Let us note that in the domain of its applicability, i. e., at An <C 1, this 
dependence means the exponential suppression of the cross section, at least 
if f(e) does not become infinite. But as An increases, the function F tree (An, e) 
becomes positive and this suppression disappears. Thus, in the latter case 
one has to take into account loop corrections to F(Xn,e) which are of the 
order (An) 2 and higher ||. 

Making use of the technique developed in refs. ||, it is possible, at least 
in principle, to find the only unknown function f(e) in (||). 

Even in the simplest case of small An the calculation of the exponent 
F tree at all energies e is rather complicated (the method of the calculation is 
described in section 0). The corresponding classical solution has singularities 
on a three-dimensional surface in four-dimensional Euclidean space (in the 
case of four dimensional theory). This surface depends on e and is determined 
in the course of the calculation. The Rayleigh-Ritz variational procedure 
enables one to obtain the lower bound on F tme . Computational procedure 
and its numerical realization are described in sections |3| and |j. The purpose 
of this paper is to explore the possibilities for the actual calculation of the 
tree exponent ([5]). In section |5]our results are compared to known analytical 
results in the limiting cases. 

2 Singular solutions and tree cross sections. 
General formalism 

Let us describe the technique to obtain the exponent for tree cross sections 
in the limit (§) 0, [|, H • We consider the process of the decay of one virtual 
particle of energy E and momentum P = into n real particles in the model 
with the lagrangian density (H). Let us write the matrix element (f3\S<p\0) in 
the coherent state representation |7J, |SL Bj, in (d+l)-dimensional Minkowskian 
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space-time ((/3| is a coherent state) 

(J3\S<p\0)= lim [ DipDip t Dip f if(E,P)e B ^ )+B f^^f )+ ^ Cdd+lx , (6) 



where 



¥>i(k) = v?k(Ti) , y?/(k) = ^(T/) , 
and the boundary terms are 

~ Jd d k^\e 2 ^ T f + Jd d kV2^e^ T f^ f (-k) 



where co> k = yl + k 2 . At the tree level the integral (|GJ) is determined by the 
value of integrand taken at the saddle point. The extremum conditions for 
the exponent are the classical field equation 

«9> + ^ + A^ 3 = (7) 

and the following boundary conditions: 

Mt) ^ a**** , Mt) ^ -§Le^ + c k e-^ 4 , (8) 

V 2 ^k 

where a k and c k are arbitrary. At t — > — oo the solution <£> c (/3*, t, x) has only 
positive frequency part. Energy conservation then implies that at t — > +oo 
the solution should not contain negative frequency parts either, i. e., c k = 0, 
and the exponent in eq. (|6|) is zero. Thus, the matrix element in the coherent 
state representation has the form of the Fourier component of the saddle 
point solution, 

A E (P*) = (P\S<f\0) tiee = f dtd^ c (/M,x)e-^ +iPx • (9) 



As follows from the coherent state formalism, the tree amplitude of the pro- 
cess 1 — > n can be obtained from the matrix element (§) in the following 

way, 

d n A E ((3*) 



Ai_> n (ki, . . . ,k n ) 



dfc . . . d(3^ 



(10) 

/3*=0 

To find the n-particle cross section, let us introduce the generating func- 
tion 

m,E) = ^ DpD(3*enp j- Jd d kp k P^A E (^(3*)A E (^(3) , (11) 

where Z is the normalization factor. The total cross section is then given by 
the following formula, 



i an 



One can check this relation by differentiating the right hand side of eq. (|TT|) 
and using eq. (0) §, [K|. 

Making use of the Cauchy formula one can rewrite the expression for cr\ r ^ n 
in the following form 



tree 

z 



If-^lJ D(3D(3*exp j-i J d d kp*(3^ A E ((3*)A E ((3) . (12) 



This integral can be calculated again in the saddle point approximation, after 
taking care of zero modes corresponding to time translations and possible 
exponentially large factors in A E (/3*). To get rid of both of them, let us 
introduce the following variables, 

In terms of these variables we have 

<p c (P*,t,x) = <p(b*,t + t ,x + x ) , 
A E (P*) = A E {b*)e iEt °- iFxo . 

Here to, xo are collective coordinates and £>k are new integration variables, 
obeying some constraint that fixes translational invariance. The form of 
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this constraint will be determined later. In terms of new variables, eq. (O) 
becomes 

a\Z e n = ^j J DbDb*dx Q dx' JA E (b*)A E (b)x 

xexp iE(t - t' ) - iP(x - x(,) (13) 
-i / A fo*fo k e^k(*o-*o)- ik ( x o-^o) _ n l n £ , 

where J contains 5-function of the constraint on 6k and the corresponding 
Faddeev-Popov; the later that does not make exponential contribution and 
will not be considered in what follows. The integration over (xq + x' ) gives 
the volume factor canceling out with Z. If there are no more exponentially 
large factors in A E , we can use the saddle point approximation (the saddle 
point of the variable (x — x' ) is equal to zero because we work in the centre 
of mass frame P = 0; we will not write this variable later on): 



a tiee {E,n) oce^- , (14) 



where W^g* 1 is the extremum value of the functional 



W tree (T } 9, b k , &*) =ET-n6- e~ e J d d k 6 k 6 k e" kT (15) 



over T = i(t Q — t' Q ), 6 = ln£, 6 k and 6 k . 

Let us now determine the constraint on 6 k . It should break the transla- 
tional invariance. We have already mentioned that we need the condition on 
6 k to get rid of the exponential factors in A£(6 k ). Let us continue analytically 
the solution of eq. (0) to the Euclidean time. Then the boundary condition 
(H) and the absence of negative frequency parts, c k = 0, implies that the 
solution should decay at Imt = r — > +oo, where r is the Euclidean time. We 
are not interested in instanton effects, i. e., we do not consider classical so- 
lutions regular in Euclidean space (in \(p 4 theory there are no such solutions 
if A > 0). Then, ip c should be singular somewhere in Euclidean space-time. 
Generally, ip c is singular on some surface r = t s (x), where t s (x) < for 
solutions smooth on the real time axis. The behaviour of the integral @ 
is determined by singularities of the function ip c , i. e., it is proportional to 
exp(_E , r m + iPx m ), where r m and x m are coordinates of the singularity closest 
to the real axis (r m < 0). Thus, to get rid of exponential factors in A E {(3*) 
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we need r m — > 0, x m = 0. This means, in other words, that we require 
the singularity surface in Euclidean space-time to touchF] the plane r = at 
the point x = 0, i. e., r s (x = 0) = 0; t s (x) < at x ^ 0. This condition 
simultaneously fixes the translational invariance (in complex time), so it is 
indeed a constraint fixing zero modes. 

So, the problem of finding the tree cross sections at any E and n can be 
formulated in Euclidean space-time and consists of the following steps: 

• Find 0(<f)-symmetric solutions <p(r, x) of the Euclidean field equations 

d 2 ip — (p — \(p 3 = , (16) 

which is singular on the surface r s (x) < 0, t s (0) = and has the 
following asymptotics at r — > oo: 

<! " X Mr, x)e" jkx = ^Le"^ . (17) 



(2tt)<V 2 ^ ' ' 

Calculate its frequency components 6 k and determine W according to 
eq. (H). 

The functional W should then be extremized over variables 6 k , 6 k (or, 
what is the same, over all singularity surfaces of the described type), 
T and 6. The tree cross section of the process 1 — >• n is then given by 
the formula (TP 



Analytical solutions of this boundary value problem can be found only 
in special cases (they will be briefly described in section |5|). Furthermore, 
in numerical computation it is impossible to extremize the functional (|15|) 
over infinite dimensional space of singularity surfaces. One can only make 
use of Rayleigh-Ritz procedure, i. e., choose some finite dimensional sub- 
class of these surfaces and extremize the functional within this subclass. 
Let us consider this process more closely. Let the functional J d d k & k ^k etJkT 
reach its minimal value for some 6 k : Jd d k fr k & k e WkT | min = C(T) > 0. Let 
us fix some family of singularity surfaces S(T). For these surfaces we have 

-'Xet us emphasize that we require only that r s approaches zero (r s — > at the point 
x = 0). 
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jd d k 6 k 6£e WkT | s(T) = C E (T) > C(T) for all T. After inserting the saddle 
point value for 9, equal to 9 = — Inn + In C(T), we get 

= n In n - n + ET - n In C(T) has an extremum at T u 
Wv(T) = n In n - n + ET - n In C s (T) has an extremum at T 2 . 

Comparing W(Ti) and Ws^) one can obtain the following inequalities: 

W(Ti) > W(T 2 ) > Wx(T 2 ) if W has a maximum at T\; 
W S (T 2 ) < W^(7\) < W(7i) if W E has a minimum at T 2 

(in the real computation the second case is realized). Thus if we limit our- 
selves to some subclass of singularity surfaces, we get a lower bound on the 
exact value of W tTe e(E,n). 

3 Expansion in spherical modes 

The following calculations will be done in (3 + l)-dimensional space-time. 
We will consider only compact singularity surfaces. 

The only requirement imposed on the singularity surface is that it touches 
the plane r = at the point x = and r s (x) < at all other spatial points. 
So, we can describe the singularity surface using the following method. Let 
us choose a sphere of the radius R s with the centre at the origin of the 
coordinate system. Let the field configuration tp be infinite at the point 
r = R s , x = and finite at all points Vx 2 + r 2 > R s . Then the singularity 
surface for this field touches the plane r = R s at x = and is contained 
inside the chosen sphere, i. e. at r s (x) < R s . This description is suitable for 
singularity surfaces of the form of a sphere that is slightly squeezed along 
the horizontal direction. These configurations are ones of primary interest, 
as we will see from the results of the calculations. We should only make 
the substitution r — > r + R s , to move the singularity to the origin. This is 
equivalent to the following change in frequency components of the field, 

K = K^ kRa , (is) 

where 6 k are Fourier components of the field singular at the point (R s , 0). 

As far as the field configuration is 0(3) symmetric, the field is a func- 
tion of two variables, <p(p, 9), where 9 is the angle between the radius- vector 
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and t axis, and p is the length of the radius-vector (in 4-dimensional Eu- 
clidean space). The Euclidean field equations can be obtained by varying 
the following action: 



1 fd(p 



1 (Dip 



2\dp) + 2p 2 \36 



S = 4tt J dej dpp 3 sin 2 9 

PminW 

Let us make use of the expansion in spherical modes, 



(19) 



p(p,9) = J2Mp)C£\cost 



(20) 



n=0 



where Cn\cos8) 



sin(n+l) 



are Gegenbauer polynomials. Asymptotically, 



as p — > oo, the functions p n (p) are of the form 

K n +i(p) 



tfn(p) = dr. 



P 



(21) 



where K n are modified Bessel functions. The coefficients 6 k of the expansion 
of this field configuration in plain waves are 



6 k = ^anV^ 



n=0 



so the integral in eq. (15) can be expressed in terms of the coefficients a n as 
follows, 



I(z) = / 6 k 6 k e" Wk ^ 3 k = 2tt 2 



n,m=0 



K 



n+m+2 



'z) - K n 



(22) 



Upon substituting the mode expansion ( |20D into eq. (fl9|) , we get the expres- 
sion for the action in terms of the spherical modes. Its extremization yields 
the equation ([16]) in terms of radial functions <p n (p)- 

The radial functions should have the form fl2"i"[) at p — >• oo, i. e., they 
should not have growing components. 

One also has to impose the second boundary condition which will ensure 
that the field becomes infinite on some singularity surface, subject to all 
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requirements mentioned in the beginning of this section. To formulate it 
precisely, one has to move a bit away from the singularity, i. e., the condition 
f(R s ,0) = oo should be substituted by the condition ip(R, 0) = A, where 
A>1/ vA- In this case one can neglect the mass term in the field equation 
near the point (-R s ,0) and approximate the singularity surface by a plane. 
Then (p in this region is 

where l(x) is the distance from the point x to the singularity surface. It is 
straightforward to see that the singularity is placed at the following distance 
from the origin, 



R s — R — 




Thus, the singularity surface satisfying the required constraints (its form 
will be described in more detail later) is determined by a set of spherical 
components c n = (p n (R), which should satisfy the condition 

oo 

J2c n (n + 1)=A, (24) 

n=0 

i. e., tp(R, 9 = 0) = A; and also the condition 

(p(R, 6 7^0) < A , (25) 

which in the simplest case of two non-zero components c n is reduced to the 
requirement that both are non-negative. 

The simplest configurations are 0(4) symmetric. They are defined in 
such a way that Co = A and c n = for all other n, and are characterized by 
only one parameter, the radius of the singularity surface R. 

Let us extremize over the parameters T, 9 and singularity surfaces. Mak- 
ing use of eqs. ( Jl8| ) and (p2|), one writes the expression ( Jl5| ) in the following 
form, 

W tree (T, 9, 6 k , b* k ) = ET-n6- I(z)e~ e , 

where z = 2R S — T. The stationarity conditions for this expression over T 
and 6 can be easily obtained, 

n = I(z)e- e , E = -I'(z)e' 9 , (26) 
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where I f (z) is the derivative of the expression ([22]). Finally, 

At? 

W tree (b k ,bi) = n \n—-n + nf(e), (27) 

lb 

f{e) = (e+ 1)T + Inl6 -\n\I(z) , (28) 
where T should be expressed through e by solving the following equation, 



which is a consequence of eq. (p26|) . 

For calculational reasons it is more convenient to perform the extrem- 
ization in a slightly different order: first fix some value of T, then find the 
minimum of I(z) over all singularity surfaces (6 k and &£) and, finally, obtain 
the corresponding value of e from eq. (|29|) . 

Let us note that the function f(e) is independent of A. Indeed, let us 
make the substitution 

Then eq. fll6|) becomes the equation for </? with A = 1, and the integral (p2]) 
transforms into 

The dependence on A in eq. (pgj) disappears. Thus one can set A equal to 
one for the calculation of f(e). 

In the case when c n c (or, equivalently, c n A) for all n > 0, it is 
possible to determine the deviation of the singularity surface of the corre- 
sponding field configuration from sphere. In this case the field (p is large at 
p = R for all 9, so in this region we can use the approximation of massless 
field. We will also assume that the radius R is large enough to consider the 
singularity surface flat at all points. Then, by making use of eq. (p3|), we 
immediately get 

AR S (9) 2 ( ' ' 



A \(p(R,9) <p(R,0). 

where AR S (9) = R s {0) — R s {9) characterizes the deviation of the singularity 
surface from sphere. The shape of a typical singularity surface found from 
our calculations is shown in fig. |1|. 
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Figure 1: Singularity surface at energy e — 10 (solid line). The 4th harmonic 
is turned on at infinity A sphere of the radius R, on which one can determine 
c n , is also shown. 

4 Numerical calculation of tree cross sections 

In the case of 0(4)-symmetric solutions (one-parameter family of singularity 
surfaces) the problem is simple: it is reduced to one ordinary differential 
equation for <po(p). It is even not necessary to solve the boundary value 
problem, one can merely fix different values of ao (i. e., initial conditions at 
infinite p) and find the corresponding singularity radii R s . The lower bound 
on f(e) obtained in this way is shown in fig. ^| by a solid line. 

If we do not limit ourselves to spherically symmetric modes only, we 
should, in principle, solve the boundary value problem (|T6D. It is difficult 
to solve it directly because one has to find (p n that rapidly decrease with 
the harmonic number n to have the sum in eq. fl2"2p convergent (and this is 
required at large p, where the field itself is small). Furthermore, in direct 
approach the configuration is determined by the shape of the singularity 
surface, where the values of the field are large, so the problem is extremely 
unstable. For the same reason one is unable to solve the boundary value 
problem formulated in terms of spherical components by fixing ip n at the 
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Figure 2: Bounds on f(e). 1) 0(4)-symmetric bound. 2) Analytical es- 
timate for small energies e. 3) Analytical lower bound for high energies. 
4) Voloshin's lower bound. 5) The bound obtained with the 4th harmonic 
turned on at infinity. 



radius R near the singularity. 

However, one can invent a method that is similar to one used for the 
0(4)-symmetric problem. Let us fix the values of a n , i. e., spherical modes 
at infinity (that is equivalent to fixing 6 k ) and solve the system of ordinary 
differential equations with initial conditions specified in this way. At some 
finite radius R the field becomes infinite (more precisely, becomes larger 
than some prescribed large number A). If the coefficients c n = (p n (R) at 
this R obey the conditions (^) and (^), then this configuration satisfies all 
requirements mentioned above. We found that this property is satisfied in a 
sufficiently large region of parameters a n , i. e., radial functions at infinity. In 
particular, one can set all a n equal to zero for all n ^ 0,k for any fixed k. 

The major problem of this method (and presumably any method that 
makes use of the expansion in frequency components) is that the number of 
spherical modes used for computation is limited. However, the errors caused 
by this truncation are significant only at short distances from the singularity 
(in the region of strong non-linearity) and do not lead to a considerable error 
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Figure 3: Bounds on f(e) (scaled part of fig. |^). 1) 0(4)-symmetric bound. 
2) Bound obtained with the 2nd harmonic turned on at infinity. 3) 4th 
harmonic. 4) 6th harmonic. 5) Analytical lower bound for high energies. 



in the determination of R s , at least for not too high energies. 

To search for the minimum of the I(2R S — T) as a function of a n (that ef- 
fectively means the minimization over 6 k , &£) we used the multidimensional 



downhill simplex method ||11|| , that does not require any additional infor- 
mation about the function which is extremized (i. e. it does not need its 
derivatives). To find the minimum, about 50 calculations of I(2R S — T) had 
to be made for different configurations. 

The calculations were performed with turning on two non-zero spherical 
modes at infinity, a and a^. All values of k smaller than 8 were explored. 
The maximum deviation from the 0(4)-symmetric result was obtained at 
k = 4. Let us recall in this respect that we obtain the lower bound on 
f(e) by the Rayleigh-Ritz procedure, so we are interested in the maximum 
value of f{s). The results for k = 4, 2 and 6 are shown in the fig. |3|. At 
smaller energies e the difference from the spherical calculation is negligible 
while at larger e the singularity radius R s becomes small and relative error 
in its calculation grows. The 0(4)-symmetric calculations were performed 
in a wider interval of energies, which is shown in fig. |2] together with some 
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analytical results (see section ffy. 

The shape of the singularity surface corresponding to the minimal value 
of I(2R S — T) is also of some interest. The singularity surface for e = 10 and 
k = 4 is shown in fig. [I]. 



5 Comparison of analytical and numerical re- 
sults 

In the limiting cases of small and large energies it is possible to implement the 
procedure of the section || analytically ||, Let us present the analytical 
results here for comparison. 

In case of small energy e one may start with the solutions of the following 
form, 

( s [2 1 



A sinh(r — r s (x)) 

which satisfies eq. ( |IED up to the terms of order 0((<9 x r s ) 2 ) (the case when 
r s (x) = for all x corresponds to the creation of particles at the thresh- 
old). One can also find corrections of the orders (<9 x r s ) 2 and (<9 x r s ) 4 to this 
expression. As a result, one obtains the following estimate ||, 

^ = !( ln ^ + 1 )--ii + is' 1327 - 9fa2) + 0(£3) • (30) 

This analytical result is shown in fig. |2] by a dashed line; it coincides with 
our numerical result at e < 0.5. 

In the case of high energies one may try to neglect the mass term in the 
field equation and consider the massless y? 4 theory. In massless theory, an 
0(4)-symmetric solution — the Fubini-Lipatov instanton — is known [|l^, [14]]. 
This solution may be used to construct the solution which is singular at the 



point r = x = [15 



^ = ^X^ + (r+ S R s r-RV (31) 

where R s is the collective coordinate determining the size of the singularity 
surface. From this solution one obtains the following bound M, |TjJ, 

/(e^oo)>lnil. (32) 

7T Z 
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This bound is consistent with our numerical calculations at e > 50. 

It is of interest to compare our numerical results with other existing esti- 
mates. An alternative lower bound on f(s), following from the direct analysis 
of Feynman diagrams, is given in ref . ]16) . As is clear from fig. |2] our bound 
is considerably more stringent. 

6 Conclusion 

In this paper we described a method for the calculation of multiparticle cross 
sections in (p 4 theory at the tree level. Our numerical results obtained in 
the 0(4)-symmetric case coincide with known limiting cases in the domains 
of validity of the latter. More precise results obtained with larger family of 
singularity surfaces indicate that even the simplest 0(4)-symmetric approx- 
imation gives nearly an exact answer. The lower bounds obtained in this 
paper are better than previous bounds. 

This work is supported in part by Russian Foundation for Basic Research, 
grant No. 96-02-1 7449a, and Soros Students Programme. The author is 
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